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We report an observation of the decay B'^ — > -DJtt^ in pp collisions at y/s = 1.96 TeV using 
115 pb~^ of data collected by the CDF II detector at the Fermilab Tevatron. We observe 83 ± 11 
Bg — » D'^Tv'^ candidates, representing a large increase in statistics over previous measurements and 
the first observation of this decay at a pp collider. We present the first measurement of the relative 
branching fraction D-tt+)/B{B° D-tt+) = 1.32 ± 0.18(stat.) ± 0.38(syst.). We also 

measure B{B+ D°Tr+)/B{B° D"7r+) = 1.97 ± 0.10(stat.) ± 0.21(syst.), which is consistent 
with previous measurements. 

PACS numbers: 13.25.Hw 14.40.Nd 



B^-B1 oscillation is expected to occur in the Standard 
Model of particle physics. Measurement of the oscilla- 
tion frequency, when combined with that for B^-B^ os- 
cillation, tests the unitarity of the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix. A deviation from 
unitarity could arise from a variety of new physics effects 



B^ meson oscillations have yet to be directly ob- 
served, with a lower limit (95% C.L.) on the oscillation 
frequency currently at 14.5 ps~^ 0. For comparison, the 
average B meson lifetime is of the order of 1 ps, and the 
B^ oscillation frequency is 0.5 ps~^. This lower limit im- 
plies that excellent proper time resolution is required to 
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Tree Color-Suppressed Exchange 



FIG. 1: Different Feynman diagrams which contribute to B — > 
Dtt decays. 

observe the oscillation. In semileptonic decays, the 
proper time measurement is degraded due to the unde- 
tected neutrino. Fully reconstructed hadronic i3° decays 
such as i?j D~Tr~^ do not suffer from this problem [Sj. 
Obtaining a large sample of such decays is an important 
first step toward measuring i?^ mixing. 

In addition, the measurement of the branching ratio 
of D-TT+ decay along with that of B° D-tt+ 

and 5+ — > D°7r+ can be used to study B meson decay 
mechanisms. As shown in Fig. H ^ D~Tr'^ decay 
proceeds only at tree level while the and modes 
have additional, non-tree, contributions. Therefore, mea- 
surements of ratios of branching fractions of these decays 
can, in principlCjisolate contributions from the different 
decay diagrams [J| . 

To date, a few B'^ DJtt^ events have been ob- 
served [alfij in e+e^ collisions at LEP. B factories, cur- 
rently running at the T(4S') resonance, do not produce 
B^ mesons. However, large samples of B^ are produced 
at the Tevatron. The ability to trigger on displaced ver- 
tices allows the upgraded Collider Detector at Fermilab 
(CDF II) to collect large samples of fully reconstructed 
5° decay modes. 

In this letter, we present an observation of 5° 
Z?77r+ decays and a measurement of the ratios of branch- 
ing fractions of B^ — > DJtt"*" and B^ D^tt"*" relative 
to the branching fraction of B'^ — > D^tt"*" decays. We 
use a sample of fully reconstructed B Dtt decays cor- 
responding to 115 pb~^ ofpp collisions at -^s = 1.96 TeV 
collected by the CDF II detector at the Fermilab Teva- 
tron between February 2002 and January 2003. Charge 
conjugate modes are implied throughout this paper. 

The components of the CDF II detector relevant to 
this analysis are described briefly below; a more com- 
plete description can be found elsewhere 0- We use 
tracks reconstructed by both the Central Outer Tracker 
(COT) and the silicon microstrip detector (SVX II) in the 
range l^y] < 1 8], where rj is the pseudorapidity defined 
as 77 = — lntan(0/2) and 9 is the polar angle with respect 
to the proton beam direction. The SVX II detector con- 
sists of double-sided silicon strip sensors arranged in five 
cylindrical shells with radii between 2.5 and 10.6 cm 0. 
Surrounding the SVX II is the COT an open-cell 
drift chamber that has an inner (outer) radius of 40 (137) 
cm. The COT has 96 layers, organized in 8 superlay- 



ers with alternating superlayers of axial and ±2° stereo 
readout. The B decays used in this analysis are selected 
with a three-level trigger system. At Level 1, charged 
tracks are reconstructed in the COT axial superlayers by 
a hardware processor, the eXtremely Fast Tracker (XFT) 
[ll| . This trigger requires two oppositely charged tracks 
with transverse momenta pt ^ 2 GeV/ c and scalar sum 
PTi + PT2 > 5.5 GeV/c. At Level 2, the Silicon Ver- 
tex Trigger (SVT) ^ associates SVX II r — ip position 
measurements with XFT tracks. This provides a precise 
measurement of the track impact parameter, do, which 
is defined as the projection of a track's distance of clos- 
est approach to the beam line onto the transverse plane. 
The resolution of the impact parameter measurement is 
50 /im, which includes approximately 30 /im contribu- 
tion from the transverse beam size. Hadronic decays of 
heavy flavor particles are selected by requiring two tracks 
with 120 /im < do < 1000 /im. The two trigger tracks 
must have an opening angle in the transverse plane satis- 
fying 2° < \A(p\ < 90° and must satisfy the requirement 
Lxy > 200 /im, where L^y is deflned as the distance in 
the transverse plane from the beam line to the two-track 
vertex projected onto the two-track momentum vector. 
A complete event reconstruction is performed at Level 3, 
and the Level 1 and Level 2 requirements are confirmed. 
Since events satisfying the displaced vertex trigger are 
dominated by the decays of promptly produced charm, 
the largest background to _B Dtt decays results from 
combining a prompt D meson with an unrelated track 
from the event. 

Reconstruction of B meson candidates begins by se- 
lecting D meson candidates. We reconstruct the fol- 
lowing D meson decay modes: Z?" K^tt^ , D^ 
K~^TT~TT~ and Dj — > 07r~ followed by — > K^K~. 
We exploit the narrow cj) — > K'^K' resonance in the 
-D7 decays to greatly suppress background by requiring 
1010 MeV/c^ < m{K+K-) < 1028 MeV/c^. No par- 
ticle identification is used in this analysis. All particle 
hypotheses consistent with the candidate decay structure 
are attempted. The track combinations which comprise a 
D meson candidate are required to originate from a com- 
mon vertex. An additional track with pT > 1.6 GeV/c 
and AR{D,ttb) < 1-5 is added and assigned a pion 
mass to reconstruct the B meson candidate. We de- 
fine AR = y^Att^'^kA^, where A(p and A77 are the 
azimuthal angle and pseudorapidity of the pion with re- 
spect to the direction of the D meson candidate. The 
selection requirements are optimized to yield the largest 
S/\/S + B for each of the decays. In the optimization 
procedure, the number of signal events (S) is estimated 
using a Monte Carlo simulation, while the number of 
background events {B) is estimated using sidebands of 
the B meson mass spectra reconstructed in data. 

We require the B meson candidate tracks to be consis- 
tent with the following constraints: the D meson tracks 
originate from a common vertex, the momentum of the 
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D meson points back, in three dimensions, to the remain- 
ing B meson candidate track, and the invariant mass of 
the D meson decay products is consistent with the world 
average of the corresponding D meson mass 01 ■ We also 
require that at least two of the B meson daughter tracks 
are consistent with the trigger requirements. The com- 
binatorial background is strongly reduced by requiring 
that L^y > 400 fim. and Lxy{B — > D) > — 150 /im, where 
the latter refers to the L^y of the D meson decay vertex 



with respect to the B meson decay vertex. The impact 
parameter of the B meson momentum with respect to 
the beam axis is required to be less than 80 /zm to assure 
that the B meson candidate originates from the primary 
vertex. The invariant mass distributions of i?^, B~^ , and 
B^ meson candidates are shown in Fig. 12 The prominent 
peak at each expected B meson mass establishes B Dtt 
decay signals, including our observation of i?° Djn^. 
We measure the following ratio: 



BiBl^^D^ ^ N{B^ eiB^ 
B{BO D-TT+) N{BO) e{B2) 



fd 
fs 



B{D- K+TT-TT-) 



B{Ds (/)7r-) • B{(t) -> K+K-) 



r 



where N{B'^) and iV(i?°) are the signal yields, fd/fs 
is the ratio of fragmentation fractions for B^ and B^ 
mesons, and e{B^)/t{B^) is the ratio of trigger and re- 
construction efhciencies. Equation ^ also applies to the 
_B+ D'^tt^ mode, with terms for B^ decays replaced 
by those relevant to 5+ decays. There are three com- 
ponents to the ratio of branching fractions measurement: 
the ratio of B meson yields obtained from fits to the 
invariant mass spectra, the ratio of signal efficiencies ob- 
tained from Monte Carlo simulation, and the product of 
previously measured production and branching fractions. 

In our analysis, we use Monte Carlo simulation to de- 
termine shapes of mass spectra and relative efficiencies. 
The Monte Carlo generation proceeds as follows. Trans- 
verse momentum and rapidity distributions of single b 
quarks are generated based on calculations using NLO 
perturbative QCD 13J . B meson kinematic distributions 
are obtained by simulating Peterson fragmentation 0| 
on quark-level distributions. Additional fragmentation 
particles, correlated 6 — 6 production and the underly- 
ing event structure are not simulated. B meson decays 
are simulated using EvtGen . The simulation of the 
CDF II detector and trigger is based upon a GEANT de- 
scription |3| that includes effects of time variation of the 
beam position and hardware configuration of the SVX II 
and SVT. 

The yield of B mesons is extracted from the invariant 
mass spectra using a binned fit, as shown in Fig. |2] 
There are three contributions to the background shape 
used in the fit. The combinatorial background compo- 
nent is modeled with a linear combination of a linear 
and exponential distribution. The shape of background 
due to Cabibbo-suppressed B DK decay is modeled 
by a Gaussian and is shown as shaded distributions in 
Fig.|2a)-(c). Backgrounds due to partially reconstructed 
B decays contribute to the low mass side of the signal 
peak and are modeled by inclusive B DX Monte Carlo 
simulation. The structures in the region close to the sig- 
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FIG. 2: Invariant mass spectra for (a) D~tv^ (b) D~tv^ 
and (c) D'^tv'^ for data, with fits overlaid. The main back- 
ground under the signal peak is combinatorial, modeled with a 
combination of a linear and exponential function. The shaded 
peak corresponds to _B — > DK decays, and is modeled with 
a Gaussian. The additional background component in the 
low mass region is due to partially reconstructed B decays. 
This background component is modeled using shapes deter- 
mined by inclusive B — > DX Monte Carlo. A sample Monte 
Carlo mass distribution for B^ D'^tt'^ events is shown in 
(d). Contributions from different sources in (d) are ordered 
by peak position, from right to left. 



nal are due to i? ^ D*7t decays, where a photon or a tt'^ 
from the D* decay is not reconstructed. The decay of 
the polarized D* produces the double-peaked structure. 
As an example of the various contributing backgrounds, 
Fig.|2Id) shows the invariant mass spectrum for B DX 
Monte Carlo reconstructed as 5+ — > Z?'^7r+. In the fit 
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to the data mass spectrum, the fractions of combinato- 
rial background and partially reconstructed -B's are al- 
lowed to float in the fit. The ratio of Cabibbo-suppressed 
B DK background to the corresponding signal is fixed 
to the world average ratio of branching fractions, with the 
trigger and reconstruction efficiencies of the two modes 
taken into account. The fitted yields of — > D^tt^, 
B+ W and B° DJtt+ decays are 1118 ± 43(stat.), 
1260 ± 42(stat.) and 83 ± ll(stat.) events, respectively. 

Trigger and reconstruction efficiencies for B"*", and 
B^ decays are determined using Monte Carlo simulation. 
The trigger efficiency differs among the three B meson 
decay modes due to differences in decay kinematics (e.g. 
opening angle distributions) which arise from the differ- 
ent masses and spins of the intermediate and final state 
particles. The ratios of efficiencies are e{B^)/e{B'^) = 
0.708±0.010 and eiB°)/e{B°) = 0.903±0.012, where the 
uncertainties are due to the limited statistics of Monte 
Carlo samples. 

World average values are used for the various 
branching fractions in Equation ^ In terms of B me- 
son production, we assume fd = fu, consistent with pre- 
vious measurement The B°/B° production frac- 
tion used in this analysis is the world average value, 
f^/f^ = 0.270 ±0.029, currently dominated by LEP mea- 
surements. The ratio of fragmentation fractions may be 
different in a hadron collider environment than in e~^e~ 
collisons. However, the ratio of fragmentation fractions 
measured by CDF [Tg] is consistent with the LEP results. 

Many systematic uncertainties cancel in the measure- 
ment of ratios of branching fractions due to the similarity 
of final state kinematics. Systematic uncertainties come 
from three main sources: fitting the invariant mass dis- 
tributions to obtain signal yields, determination of the 
ratio of efficiencies, and uncertainties on external inputs. 
The contribution to the systematic uncertainty from ex- 
ternally measured quantities is calculated by propagating 
world average uncertainties. Systematic uncertainties on 
the signal yields are determined by comparing the fitted 
yields after changing the invariant mass region of the fit 
and varying the background shape within the range al- 
lowed by the B — > DX Monte Carlo statistics and world 
average uncertainties on the branching fractions of par- 
ticipating B decays. Systematic uncertainties on the ra- 
tio of efficiencies come from physics sources such as the 
choice of pt spectrum or meson lifetime, and detector 
sources such as inaccuracies in the XFT and SVT hard- 
ware simulation. The uncertainty due to a given source 
is estimated by the shift of the ratio of efficiencies when 
the effect of that source is modified in the Monte Carlo 
simulation. The effect of the choice of B meson pT spec- 
trum is estimated by re-weighting the Monte Carlo to 
match the measured B hadron pr spectrum 0. Since 
the trigger and analysis selection only accept events in 
which the B meson is displaced from the primary in- 
teraction point, the efficiencies depend upon the B and 



D meson lifetimes. To estimate uncertainties due to B 
and D meson lifetimes, the Monte Carlo is re-weighted 
with different lifetimes within the world average uncer- 
tainties. Due to the different specific ionization of tt^ and 
in the COT, kaons are 6% less efficient in satisfying 
the XFT requirements ,18|. The Monte Carlo simula- 
tion is re-weighted to reproduce this effect. To estimate 
the uncertainty due to imperfections in simulating the 
signal selection requirement efficiencies, we compare ef- 
ficiencies between Monte Carlo simulation and sideband 
subtracted B~^ and B'^ signal for every selection require- 
ment separately. In the case of the B^ / B^ branching 
fraction ratio measurement, there is an additional sys- 
tematic uncertainty associated with the fact that the B^ 
decay has four tracks in the final state while i?+ has 
only three. We infer the corrections and systematic un- 
certainties due to the fourth track by comparing trigger 
and reconstruction efficiencies in data and Monte Carlo 
for semileptonic B — > D*^l+X, D*~ l3°7r~ decays be- 



tween the K^TT and A'+7 



final states. 



The systematic uncertainties are summarized in Ta- 
ble ^ The total systematic uncertainty is obtained by 
adding individual contributions in quadrature. Using 
Equation^ we obtain the following values for the ratios 
of branching fractions: 



g(i?g ^ D-7r+) 
B{B° D-TT+) 



B{B+ D°Ti+) 
B{B" D-7r+) 



= 1.32±0.18(stat.)±0.10(syst.) 

±0.34(BR) ±0.14(PR) (1) 



= 1.97±0.10(stat.)±0.15(syst.) 

±0.14(BR) (2) 



where BR and PR refer to the uncertainty on the ra- 
tio of D meson branching fractions and i?^ meson pro- 
duction relative to B^, respectively. Under the assump- 
tion of isospin invariance, our measurement of the ratio 
B{B+ D°'!T+)/B{B° D-Tr+) is consistent with the 
world average j2i]. This provides a high statistics cross- 
check of the measurement procedure. 

In conclusion, we have presented the first observation 
of 5° D~7T^ decays in pp colhsions, and the first mea- 
surement of the Bg D^TT'^ branching fraction relative 
to the D~Tr'^ branching fraction. The precision of 

this measurement is currently not adequate to separate 
the contributions of different decay diagrams We 
expect the measurement precision to improve as world 
average values of D meson branching fractions and B 
meson production fractions improve. 
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